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ABSTRACT 

The most luminous X-ray source in the Local Group is associated with the nucleus of M33. This 
source, M33 X-8, appears modulated by ~ 20% over a ~ 106 day period, making it unlikely that the 
comb ined emission from unresolved sources cou ld explain the otherwise persistent ^10 39 erg s _1 X-ray 
flux (Dubus et al. 1997, Hcrnquist ct al. 1991). We present here high resolution UV imaging of the 



1997| , pcrnquist ct al. 199l| ) 
nucleus with the Planetary Camera of the Hubble Space Telescope undertaken in order to search for the 
counterpart to X-8. The nucleus is bluer and more compact than at longer wavelength images but it 
is still extended with half of its 3T0 38 erg s _1 UV luminosity coming from the inner 0'.'14. We cannot 
distinguish between a concentrated blue population and emission from a single object. 

Subject headings: galaxies: individual (M33) - galaxies: nuclei - Local Group - ultraviolet: stars 



1. INTRODUCTION 



The nearby galaxy M33 hosts the most luminous steady 
X-ray source in the Local Group, X- 8. This source, with 



Lx ~10 39 erg s 1 ( |Long et al. 1996 ), is coincident to 



ithin 5" with the nucleus of the galaxy ( [gchulman fc| 



Brogman 1995 }-. — Difforont mo dole woro invoked for X 8 



including a quie scent mini-AGN ( Trinchicri et al. 1988| 
Peres ct al. 1989] ) , a collection of X-ray binari es (|Heruquist 



et al. 1991|) ~and a new type of X- ray binary (pottwakl, ct] 



al. 1987| ). Our ROSAT studies ( pubus et al. 1997| ) have 



shown that X-8 is very steady on both short and long time 
scales, except for low amplitude (~20%) variations which 
appear modulated on a ~106 day period. This strongly 
favors a single source explanation for X-8. 

We have interpreted the modulation in X-8 as "superor- 
bital" , similar to that seen in a number of bright galactic 
X-ray bi naries which were monit ored by e.g. the Vela 5B 
satellite (Smale & Lochner 1992). X-8 is then likely to be 
a >10Af© black-hole X-ray binary (the high mass is re- 
quired to account for the observed luminosity) but with 
a companion in a much shorter orbital period than 106d. 
This is supported by the extremely low velocity dispersion 
of the nucleus which limits the mass of a central black 
hole in M33 to <5xl0 4 M Q ( [Kormendy fc McClure 1993 , 
KM93). This and the high central stellar density imply 
that the nucleus is an extremely relaxed, post core-collapse 
stellar system (comparable for instance to a galactic glob- 
ular cluster such as M15). A significant number of stellar 
collisions/interactions could have taken place, eventually 
leading to the c reation of exotic interacting binaries (e.g. 
Hut et al. 1992j ). 



The next step toward unravelling the mystery of X-8 
would be to identify its optical counterpart. However, even 
with optimistic Lx /L opt ratios for either X-ray binaries or 
AGN the counterpart would only have V~21 compared to 
a core brightness of V~14. But with the optical spectral 



type of an F supergiant, the dominance of the M33 visual 
core cannot extend to UV wavelengths where the hot / flat 
spectrum of X-8's associated disc ought to be a significant 
contributor. This is true des pite evidence for a color gra- 
dient in th e nucleus (KM93, [Mighell fc Rich 1995|, |Lauei 



et al. 1998) suggesting that a period of recent star forma- 



tion has taken place and/or that collisions have modified 
the central star population. Here we report an attempt to 
find the counterpart using the UV imaging capabilities of 
HST. 

2. OBSERVATIONS AND REDUCTION 

Observations were carried out with the HST WFPC-2 
on June 12, 1997 using three different filters. The nucleus 
was positioned at the center of the Planetary Camera (a 
1:33:51.1, S 30:39:39, J2000). During the first orbit two 
1200s exposures were made with the F160BW ('UV filter, 
A1491). In the following orbit, two 800s exposures were 
made with the F300W filter ('U' filter, A2942A)_and one 
500s exposure with the F439W filter ('B' filter, A4300A). 
All these exposures were made with the gain setting at 7. 

In addition, we have extracted recalibrated archival data 
from the Space Telescope European Coordinating Facil- 
ity (ST-ECF) Archive. These data included two 40s ex- 
posures with the F555W at gain 14 ('V filter, A5397A), 
two 40s exposures with the F814W at gain 14 (T filter, 
A7924A) and six 300s exposures with the F1042M at gain 
7 (A10190A) filters, all centered on the nucleus and dating 
from September 26-27, 1994 . The archival V and I data 
were previously discussed by Lauer et al. 1998 (L98). Fig. 
1 shows the central region of the reduced UV, U and B 
images. 

The data were reduced using the HST calibration 
pipeline. The signal-to-noise ratio of the only existing 
F160BW flat was quite low. Following advice from the 
WFPC-2 group at STScI, we decided to use the F255W 



2 



flat field. Effects of cosmic rays were reduced on those im- 
ages for which we had multiple exposures with the IRAF 
stsdas routine crrej. Images of the nuclear region as ob- 
served through the F160BW, F300W and F439W filter 
with the PC are shown in Fig. 1. 

3. ANALYSIS 

Our values for the total flu x of the nucleus in t he differ- 
ent bands agree with those of Gordon et al. 1999. The ra- 
dial profiles of the B, V and I data are also consistent with 
KM93 and L98. As is apparent from Fig. 1, the nucleus of 
M33 appears more concentrated in the F160BW filter than 
in the longer wavelength filters. Nevertheless, as shown in 
Fig. 2, the profile in the F160BW image of the nucleus is 
extended compared to the profile of the star located about 
1" NNW from the nucleus and to PSF profiles ca lculated 
using the H ST PSF generating routine TinyTim (Krist & 



Hook 1997| ). Further investigation of the radial structure 
of the UV emission calls for deconvolution of the data tak- 
ing into account the different instrumental effects in the 
Planetary Camera (L98). Since the UV image does not 
have a large enough signal-to-noise to allow for a proper 
deconvolution, we chose instead to fit convolved models to 
the data. 

3.1. Fits of extended emission models 

Following KM93, we fitted radial profile models of the 
form: 



E (l + (r/r ) 2 )" 



(1) 



The model and the PSF are oversampled on a 4x4 grid 
for each PC pixel. For a given set of (r ,n), the convolved 
model is moved on the 4x4 grid, rebinned to the PC resolu- 
tion and compared to the data. The comparison with the 
data is performed in a 64x64 pixel aperture (about 3"x3") 
but we have verified that larger and smaller apertures gave 
similar results. We have assumed an A type spectrum for 
the PSF but other choices do not affect our conclusions. 
The parameter r Q is varied between 0.05 and 2 PC pixels 
and n is varied between 0.5 and 2. The results from the x 2 
minimization routine are presented in Tab. 1 and Fig. 3. 
The quoted errors correspond to 10% higher values than 
at the minimum of the fit function. In the noisy UV band 
only lower bounds to the parameters could be extracted. 
Here we give errors on the FWHM instead of on r Q . 

The FWHM of the models decreases in accordance with 
the blue colour gradient. In all cases n w 1 suggesting that 
the distribution of light at large radii is the same in all the 
filters. This is consistent with the flat colour profiles ob- 
served for r 3> r a . With n fixed at 0.75 as in L98 we also 
find a FWHM for the F555W data of 0'.'07. This is clearly 
not the best solution when n varies (Tab. 1). Wc find 
a higher FWHM (0'.'09). We note that KM93 find n be- 
tween 0.8-1.3 and a FWHM below O'.'l. The similar values 
found for r Q in the V, I, F1042M and (to a lesser extent) 
B bands indicate that their radial profiles are comparable 
(i.e. the colour gradients are much reduced between those 
bands than when compared to the UV and U so that, for 
example, on first approximation the V-I colour gradient is 
negligible when compared to UV- V) . 



3.2. Fits with an additional point source 

Since the V, I and F1042M show very close light dis- 
tributions, we investigated whether the compact emission 
from the UV and U filter could be explained by an under- 
lying extended population having the V band distribution 
plus a blue point source. We fixed r Q =l PC pixel, n=l 
and superposed at the center of this model a point source 
of varying relative strength. For the V and redward filters, 
the best fits were consistently obtained with a nil contri- 
bution from the point source. However, a point source 
of increasing strength was needed in B, U and UV. Only 
those models with fit values within the errors of the previ- 
ous extended emission fits (Tab. 1) were kept i.e. the fits 
here are as good or better than the previous ones (see Fig. 
3). The contribution of the point source to the total flux 
within a l'/45 circular aperture is summarized in Tab. 2. 
The corresponding magnitudes are given in the VegaMAG 
system usin g updated tables for the zeropoints (Holtzman 



et al. 1995). 



The best fits for the B, U and UV are shown in Fig. 3 
by dashed lines. They are indistinguishable from the ex- 
tended emission fits. As a result we cannot, based on the 
data in hand, distinguish between a model in which emis- 
sion is extended in the UV but has a smaller core radius 
than at longer wavelengths and a composite model con- 
sisting of a point source and an underlying distribution 
characterised by the visible light profile. 

4. DISCUSSION 

The nucleus of M33 has a composite spectru m ranging 
from A7V at A-3800A to F5V at A-4300A flO'Comici: 
1983). This requires at least a two component popula- 



tion in most models with the blue emission being due to 
young metal-rich stars. The colour gradient in B-R im- 
plies this young population is more centrally condensed. 
The nucleus could have been the site of episodic starbursts 
with t he youngest stars being about 10 Myr old ( O'Comicl"] 
1983| , |van den Bcrgh 199ll Schmidt ct al. 1990| and refer- 



ences therein). Recently Gordon et al. 199E have argued 
that a single 70 Myr old starburst reproduces the UV to 
IR spectral energy distribution within 4'.'5 of the center 
if dust is correctly taken into account 1 . The nucleus is 
also very similar to a glob ular cluster and is likel y to have 
undergone core-collapse (Hcrnquist et al. 1991, KM93). 
Blue stars formed in collisions at the center might explain 
the colour gradient. This model may have difficulties ac- 



1991 



coun ting for the UV lum inosity (Hernquist et al. 
L98, ICordon et al. 1999| ). 
1 Massey et al. 1996" detected the nucleus in the UV but 
with 5" resolution. Hence they had proposed that the 
blue component of the nucleus could be due to unresolved 
emission from a few hot stars. However, the HST data 
shows the UV emission is very compact with no stars of 
comparable brightness within 5" of the center. The total 
flux is about 6.3T0 -15 ergs-s _1 cm~ 2 A _1 in a 4'.'55 aper- 
ture or about 2.8T0 38 ergs-s _1 at 800 kpc. From the best 
fit model, ^50% of the UV light comes from the inner 
0'.'14 of the nucleus. Any model of the structure of the nu- 
cleus has to explain this UV emission from a region only 
~ 0.55 pc across. If a point source is present, this source 



e m^i 



ordon ct al. 1999 propose that X-8 is a high mass X-ray binary with an early B companion. But as had been noted by jO'Connell 1983 



the high mass tranter rate needed to power the 10 39 ergs s 1 luminosity implies an uncomfortably short evolutionary time-scale (~ W'° years) 



fern 



3 



is responsible for —30% of the UV flux within 0'.'14. The 
contribution from the underlying extended population is 
subsequently reduced. 

The nucleus of NGC 205, at a comparable dis- 
tance of 720kpc, is in many ways similar to that of 
M33, but without the X-ray source. It is globu- 
lar cl uster like, has a comparable My and a low up- 



per Unfit of 9-10 4 M,r ) on the mass within the cen- 



tral pc ( Heath Jones et al. 1996| ) 



These authors 

find that the nucleus is more extended, with a F555W 
FWHM of 0"2, ^ F55 5w=1.8-10- 15 ergs-s^cm^A" 1 
and only Ffi6obw=6.0-10 -16 ergs-s _1 cm _2 A _1 (within 
0'.'273). Using the same aperture on the M33 data we 
find for M33 F F555W =3.4-10~ 15 ergs-s^cm^A" 1 and 
-FpifinRw=2.5-10~ 15 ergs-s _1 cm _2 A _1 . Excluding the 



F555W = -' i " 1 >' 

Fi60BW=2.5-10" 15 ergs-s _1 c 
contribution from the point source, the F160BW flux of 
M33 is about 1.5T0 -15 ergs-s _1 cm _2 A _1 and the ratios 
of the fluxes between the bands become comparable. 

If there is in fact a point source at the center of 
the nucleus of M33, the magnitude estimates in Table 2 
are consistent with a Rayleigh- Jeans tail. We estimate 
the total UV flux of the source would be — 1.0T0~ 15 
ergs-s _1 cm _2 A _1 , which would suggest L opt / 
which is reasonable for the type of X-ray source we have 
postulated X-8 to be. This source would be responsible 
for most of the colour gradient in UV-B and U-B and thus 
the very compact appearance of the nucleus at these wave- 
lengths. As it contributes only -18% of the total F160BW 
flux and —8% of the F300W flux, the spectral energy dis- 
tribution of the nucleus is not changed much and the con- 
sequences on population synthesis studies should be minor. 
However, it does have a major influence in that it relaxes 
the constraints on e.g. mass segregation to explain the 
very strong colour gradients in UV-B and U-B. 

The star located about 1" from the nucleus in M33 has 
5=19.45, V=19.25, M F3 oow=18.05 and M F i 60 bw=17.60. 
The count rates are compatible with an AO type spec- 
trum which would make it simil ar, although fainter, to the 
two A supergiants detected by Massey et al. 1996 . The 
fluxes in the F300W and F160BW bands, at the distance 
of M33 (-800 kpc), are - 10 37 erg-s" 1 . The positional 
accuracy of the ROSAT HRI does not rule out this star 
as a possible counterpart to the X-ray source X-8. The 



ratio Lx /L opt — 100 and the observed U — B and B — V 
would a gree with what is expected from a low-mass X-ray 
binary ( |van Paradijs fc McClintock 1995 ). But its abso- 
lute magnitude (My «-5.2) would make it more similar 
to a high mass X-ray binary. The main argument against 
this star being the X-ray binary is X-8's unique c haracter 
and special location at the nucleus. Given that [Massey 
et al. 1996 found —300 analogous UV sources in M33, it 



would be remarkable that the one near the nucleus is the 
most luminous X-ray source in the Local Group. 

5. CONCLUSION 

The UV high resolution images obtained with the HST 
Planetary Camera show the nucleus of M33 is extremely 
compact. We have fitted convolved models to the ra- 
dial profiles in the different bands from which we find the 
FWHM of the nucleus in UV to be ~0'/035 and -0"090 in 
V. About half of the UV flux comes from the inner 0'/14 
of the nucleus. The UV and U profiles are also well fitted 
if one assumes a blue point source superposed on an ex- 
tended population with the same FWHM as in V. If this 
is the correct model for the nucleus, then this point source 
is likely to be the UV counterpart to the very luminous 
X-ray source X-8. Such a counterpart would be responsi- 
ble for most of the strong colour gradient seen in UV. Its 
contribution to the total UV flux of the nucleus would be 
about 18%. Models for the structure of the nucleus still 
need to account for the remainder of the UV flux but the 
constraints on population segregation (more compact blue 
star population) are reduced. High spatial resolution UV 
spectroscopy of the nucleus is the obvious next step, which 
we will undertake shortly. 
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Tabic 1: Fits of extended emission models 



Filter 


r Q n 


FWHM 




(PC pix.) 


(tf.'OOl) 



F160BW 


0.4 


1-1-0.2 


35_25 


F300W 


0.5 


i-o±g:? 


45±f 5 


F439W 


0.8 


i-o±g:i 


yr+20 
10 -15 


F555W 


1.0 


i n+0.05 
1 - u -0.05 


90 +1 5 5 


F814W 


0.9 


9+ 01 


90±« 


F1042M 


1.0 


9 +al 


lOOt^? 



Table 2: Fits with additional point source 



Filter 


Point 


source flux & magnitude 




within l'/45 aperture 




(%) 


(VEGAmag) 


F160BW 


18+4 


16.9 


F300W 


8t* 


17.5 


F439W 


It? 


19.7 




Fig. 1. — Central portion of the Planetary Camera showing the nucleus and the nearby star in the F160BW 'UV, F300W 'U' and F439W 
'B' filters. Each pixel is 0.0455" i.e. each image is about 2.3" X 2.3". 
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0.1 0.2 0.3 
R (arcsec) 

Fig. 2. — Radial profiles of the nucleus and nearby star in UV. The nucleus is shown by the filled circles + continuous line. The star NNW 
of the nucleus (see Fig. 1) is shown by the triangles + continuous line. The best fitting PSF to the star is shown by the dashed line. The 
nucleus in UV is clearly extended. 
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F1 60BW 
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F300W 
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0.0 ri ■ ■ ■ ■ i 



0.1 0.2 0.3 
F439W 



0.1 0.2 0.3 
F555W 



1.0 * 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




0.0 h 1 1 1 1 1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



0.1 0.2 0.3 
R (arcsec) 



0.1 0.2 0.3 
R (arcsec) 



Fig. 3. — Radial profiles of the nucleus in the different filters (data points). The best fitting PSF convolved models assuming only extended 
emission (paragraph 3.1) arc shown by straight lines. The unconvolvcd models (Eq. 1) arc shown as dotted lines. PSF convolved models 
assuming extended emission and a point source (paragraph 3.2) are shown with dashed lines in the F160BW (UV), F300W (U) and F439W 
(B) profiles. Normalisation factors are (in Mag/arcsec 2 ) 11.67 (F160BW), 11.65 (F300W), 12.10 (F439W) and 11.76 (F555W). 



